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I. Introduction 

of^thInTtTi*fif reflectance and transmittance 
w T . ^ common method of obtaining 
the optical constants n amd k of the metal, '-^ In 

Tnnr^J^^-J^^^'^T^'^^ "«"^"y been made 
I JZT '"'"'^^'"^■r Our purpose is to describe some 
genera procedures for analyzing measurements made 
at obhque mcidence. Solution of the Maxwell- 
2nnlT b.7'?dary-value problem with oblique inci- 
dence to obtain reflectance R and transmittance T in 
^ terms of the material constants n and k is straight- 
forward,' but the resulting functions are very compli- 
fo oh;„- of inverting the functions 

t diffi.?u °f measured and T 

^ difficult This difficulty has presumably . dis- 
couraged attempts to make much practical \i8e of 
comnntT '"l" ""^^"^ements, but we shall show how 
fTlf L ^'•^P'^^'ig a.nd inversion techniques make 
feasible the analysis of thin-film measurements at an^ 
inH^nf ''"''f^'''^- I"^, addition, a judicious choice of 
ipc dent angle and polarization can cope with some 
troublesome regions of n-k space where normal inci- 
dence measurements are inaccurate." Some of these 
results and conclusions have been obtained previously. ^.o 
bounnS l"'"'"/''^' T^'OP'^' homogeneous metal, 
PnT ^ ^-.uP'T^P*^^"^^ surfaces, without being 
concerned with the validity of these assumptions fof 
real evaporated films." We exclude from considera- 
tion monocrystalline films of anisotropic metals this 
limitation is not of importance at present because ex! 
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perimentally it is very difficult to deposit single crystal 
layers and because the theoretical analysis of the 
optical measurements would be inordinately compli! 
cated. Therefore the problem is, starting with eqSa- 
tions for R and T in terms of n and /under tE 

tTonTof/^and ' *° '^^"^^^ ^^^^^'^'^ '^^^^ 
tions of R and r for measurement and second to invert 

the equations to obtain values for n and k 

r, .^aT^^'^'^a'^ ^^^"^^^ * ^ depend not only on 
n and according to electromagnetic theory under the 

tt^^'^lm thickness of the 
thm film and on the wavelength, angle of incidence 
and polarization of the incident light. The £ tS 
quantities are usually determined by the disposition S 
the measuring apparatus. The values of n and Tare 
expected to depend on wavelength of the light (photon 
energy which excites the electronic transition), but not 
on polarization and angle of incidence (in the isotropic 
c^e) or on thickness, since n and k are material con- 
stants independent of the bounding surfaces. Thert 
Wth" Z^"' i° 1^ and A, a function of wavj 
length, at each wavelength one must measure thret 
values from among the foUowing possibilities: (1) R 
at specified incident angles and polarizations, (2) T 

?i fil^'1v/''"'^^"l^^?^^' ^'^'^ polarizations, and 
(3) film thickness. _ Varying the thickness is possible 
but inconvenient since separate samples must be fabri- 
cated (It also turns out to be less accurate," exceot 
in certain regions of n-k space which are insensitive to 
other parameters.) Measurement of film thickness 
can be made by means that are independent of the 
!tJ "\easurements, but the accuracy is usuallv 
not better than about ±30 A except under especiX 
careful conditions.' One may therefore consider 

r\^Tt (f^'^'I^T^. possibilities for 

R and T (from which the thickness could then be de- 
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termined as well as n and fc). We wish to choose opti- 
mum incident angle and polarization. ,j u J- 

The'limiting case of infinite thickness should be dis- 
cussed separately, since it is often used experimentally, 
but of course the thickness disappears as a parameter 
and the transmittance is always zero. One can there- 
fore make two measurements of R at different polariza- 
tions or incident angles on a bulk sample. An en- 
tirely different approach involves making a single 
measurement of R (usually at normal incidence) and 
using a Kramers-Kronig relation^ to obtain the two 
quantities n and fc. The latter approach has the ad- 
vantage that it can be extended farther into the vacuum 
ultraviolet region where polarization of the incident 
light is difficult to control. It has the disadvantage 
that one must make measurements in this region, even 
in order to determine n and k only in the visible region. 
. One must also make an assumption about the wave- 
length dependence of the reflectance outside the mea- 
surable range. Indeed, part of the interest m the 
techniques discussed here is to check the validity ot 
such assumptions. Finally we mention the Drude 
(ellipsometry) method, in which one measures the rela- 
tive amplitudes and phases shifts of the reflected 
polarized light instead of its intensity. (Bulk-sample 
techniques have been used to mea^sure anisotropic 
effects in single crystals.") 

We shall discuss the behavior of R and T under dit- 
ferent conditions of polarization, incident angle, and 
film thickness, in order to facilitate a choice of the best 
combination of measurements on thin films. The 
choice depends on the examination of contours plotted 
in n-k space. This procedure has been used previously 
for infinite thickness^** or for normal incidence.* Then 
we shall describe a computer-based iteration scheme 
that inverts the functions R and T to give the optical 
constants of the metal studied, and give an example to 
show how the scheme works. Comparisons will also 
be made with some bulk-sample methods in terms of the 
same criteria. We do not attempt a comprehensive 
evaluation of instrumental limitations of various 
methods,^ which would give an over-all assessment of 
the best^method for each circumstance. 

II. Theoretical Formulas for R and T 

A hght wave is assumed incident from a dielectric 
medium 1 (e.g., air) onto the plane surface of metalhc 
medium 2 at an angle of incidence di with the normal to 
the bounding surface. The back boundary of the. 
metal layer is a plane-parallel interface with another 
dielectric medium 3 (e.g., a glass substrate). The re- 
flectance and transmittance^^ of ^he metal layer are 
given by (tilde denotes a complex quantity) 



The Fresnel coefficients for perpendicular polarization 
of the electric vector are 

fij. = (m COS01 - ns cos5a)/(ni cos^i -h n* cosSa), 
= 2ni cosdi/(ni costfi -h n% cosSt), 
and for parallel polarization, 

f„p = (n, cos^i — ni cos52)/(na cosOi -h ni cos^O, 

(ij** = 2ni cosdi/(n2 cos^i 4- ni cos^^a). 

(The permeability of all three media is assumed equal 
to MO for free space.) 



where 



R = T = (nM\l\\ 

f = [f,j -H exp(i«)]/[l + fiifti exp(i6)], 
I = IWjj exp(i«/2)l/[l + fi2**23 exp(i5)]. 



d = 4ir(d/\)n2 cos^a, 

where d is the film thickness and X is the vacuum wave- 
length of the incident light. Similar expressions are 
obtained for and t2z by permuting the subscripts; 
713 and Bz are real. ^ t i . 

The angle Ot is the angle of incidence of the light 
wave and ni is the refractive index of medium 1, The 
complex angle 62 is defined by Snell's law, sin(92 = 
m sini9i ; also = n + ik, where n and k are the opti- 
cal constants of the metal that we wish to find. We 
make the definition 

ni cosSi = w 4- 

Then one finds^ 

= i{(n« - A;* - «i* sin^di) 

-f [4n2fc2 -I- (n» - Aj« - Wi* sm»0i)«lM. 

1,8 :^ i{ -(n* - A;* - ni* sin^^i) 

-h [4n2A;* -h (n* - A;« - sin=^i)*]*}» 

cos#2 - (w -t- w)/(n -h ik). 
Finally, 



The above formulae for u and v appear somewhat 
simpler when expressed in terms of the complex dielec- 
tric constant i = €r + iu. where e = In addition, 
€ is more directly related to the electron theory of the 
metal. Nevertheless, we present our discussion in 
terms of n. The mapping e = is a conformal map- 
ping, and so all the conclusions that we draw below 
about the intersections of contours of constant R and 
constant T in the n-k plane would hold also in the 
er-u plane. The former has the advantage of magni- 
fying the interesting region hear the origin. 

These formulas are sufficient to determine R and 1 
in terms of n and fc. In the application to metals we 
should note that the derivation from Maxwell s equa- 
tions assumes local conductivity. If nonlocal effects 
are included, a longitudinal wave, gives anomalous ab- 
sorption of p-polarized light near the plasma fre- 
quency, ^^-'^ but this is normally very small. Ihe 
ordinary transverse wave, excited by s-polarized light 
as well, is modified by the anomalous skin effect, 
which gives larger anomalies in the infrared, especially 
if the electron reflection at the metal surface is diffuse. 
Apparently, however, the above formulas can be used 
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0 < < 4 ^M^' e^ch curve n ha^ the value indicated and 
" S * < 4. (b) For each curve k haa the value indicated and 
0 < 1 < 3. 



^Wn lff ? *he optical constants, and the anomalous 
8km effect can be introduced afterwards a modific^^ 
tion m the classical Drude free-electron theory 

III. Contours of R and T 

y^Z ZV^Tf ? *° *° ^^^^^"^ experimental 

values of n and A; from measurements of R and T 
According to the formulas of Sec. II, R and T are 
single-valued functions of n and k, bit^it ?urns out 
Thtrf ^'"'^■^ "^"Itiple-valued functions of 22 andT 
Iherefore it ,s essential to gain some insight into the 
functions before attempting to invert them. Errone^ 
ous solutions have been atfributed to unawareness o^ 
this problem." To illustrate this point, in Fi? 1 
we plot contours of constant n and k in the R-T pLe, 
of andT iiicidence.- Only certain^pairs 

of R and T values are physically possible. They must 



be below the Ime -f- T = l, and above the (common) 

d/\ and 0) Two solutions for n and k are possible 
over much of the allowed area, for a particular R and T 

J.uJZf' '^fr * considered, multiple 

solutions would occur over all of the allowed area. 

plane and comparing our results with previous work 
and ?r thf ' ^ ' t° ^* contours of constant I 
nat*.Iv ? ^^^''^ °^ ^ie- 1- Unfortu- 

Tfl !5 ca'^not be inverted analyticaUy, 

andso the contour-plotting routine is more costly in 
computer time. The functions R and T were cSci- 
lated over a fine grid of points in a selected region of the 

tion, and d/\ If the calculated value of the func- 

T?S plotted (chosen equal to an assumed error of 1% 
of fuU scale), the point was plotted before we went on 
to the next grid point. For the contours shown here 
we have used n, = 1 (air), = 1.46 (quartz) ' 
Kesults for normal incidence are redrawn in Fie 2 
for d/X = ai. (The width of all the folWng Con- 
tours IS =-0.005 for and for T > 0.1, and ^0 OOlTr 
L^^r. . y ■ f '^^^^d intersection sug- 

fnd A- V r'''? ^^P«""^enfc^I error in the inferred n 

and I values. The absolute error in n and k is lar^e 

hrst feature is much more serious, since the region of 

(esp cfalfv W tt^''" °f measurement 

(especially for the n value, m this particular example). 




,1 - O l- The contour width mdicates a nominal error of 
mea..urement, to show the uncertainty in the derived « and * 

values. 
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Fie ' 3 Contours of const. R for normal incidence and d/X 
0.05, 0.2. 




Shows R contours for d/\ = 0.05 and 0.2. These may 
be useful, especially where the R and T contours are 
tangent, but we conclude (in agreement with Nilsson) 
that variation of film thickness probably do^ not give 
the most satisfactory method of analysis because of 
the uncertainties in thickness measurements. 

To give some orientation toward the region ot inter- 
est for real materials, in Fig. 4 we sketch on the same 
n-k plane values from the literature for some repre- 
sentative materials. Simple metals like Na avoid the 
difficult region of tangency, but others, notably Ag, 
enter it. Also the semiconductor Ge crosses it m the 
visible spectrum. Thus there is reason to consider 
the oblique incidence contours, even aside from their 
role in eliminating ambiguities and determming d, 
to be discussed below. 

We begin with two general observations which will 
facilitate limiting our discussion to the interesting 
cases. First, at 60° angle of incidence a plot of ft 
and T contours gives intersections that are qualita- 
tively similar to Fig. 2 for normal incidence, for both 
s and p polarization. That is, the region of tangency 
shown in Fig. 4 and the multiple intersections still 
occur— in about the same location for s polarization 
and shifted roughly 0.5 to the right for p polarization. . 
Therefore, we do not make a major gam by changing , 
the angle, for an R-T measurement Second, the 
s polarization contours for both iJ and T are almost 
exactly parallel to the respective normal-incidence > 
contours up to d/X = 0.2, even at an incident angle , 
as large as 75°. At a particular n and fc, of course, . 
the wHue of the reflectance is higher and the trans- 
mittance lower the larger the incident angle. There- . 
fore, except for considerations of convenience and 
instrumental accuracy, the «-polarization case may be 



Fig 4 Experimental values for some representative materials. 
The cross-hatched area is where the normal-incidence « and 2 
contours are tangent for values ot d/\ from 0.05 to 0.2. 

The points of tangency lie near the line fc = n, as has 
been observed previously.* Another related difficulty 
stems from the multiple intersections of a pair ot 
contours. This ambiguity cannot be resolved without 
additional information.' . j _ n 

Normal incidence contours plotted tor d/\ - u.uo 
and 0.2 show the same general features as tig. ^■ 
The tangency of « ami T contours still occurs near fc - 
n, and multiple intersections occur. All the 1 con- 
tours for different d intersect at relatively small angles, 
as previously remarked," giving large uncertainty m n. 
The R contours, however, intersect at angles greater 
than about 30° almost everywhere, and multiple inter- 
sections do not occur in the range considered. * igure 6 




Fig. 5. Contours of const. R for « and p polarization with 60° 
incident angle and d/X = 0.1. 
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Fig. 6. Contours of const. T for , and p polarization with 60- 
meident angle and d/\ = 0.1. 

equated to the normal-incidence case in our diacussion. 

ihe p-polanzation contours, however, differ sie- 
nificantly xn shape from the s-polarizatior^ (or normll 
mcidence) contours. At 60» incident angle th^se 
pairs of i2 contours (Fig. 5) and T contourf ^Fig 6) 
intersect at useful togles in the region n < 2 k < 2 
but outeide this region they are approximately parallet 

I?8o the r^^"^ ^^S^^'^ extends further 

J. I, .f ^1^°'' ■^^^S^'' smaUer d/\ is. (We 

shall allude to the latter point in connection with our 

TwiT'^'' measurements on bulk samples ) 
Sn 15^! measurement of i? or T at large angles on 

rS^ Tb 5?^?"°" "^^'"^ ^^d'*'^'^^ information in 
me range 1-6 eV for many materials. 

. Jixperimentally, when a spectrophotometer is used 
e'a^rThSr ''-^ various . angles arrLpecTaLy 
aw' 7 Iv T ^'""^ ^^^'"Ple" for the case of Ag 
' «\vV measurement of T at 60° with 

^polarization, m conjunction with normal-incidence 
feJweeT ^'^^^t^' ^^^only resolves the ambiguity 

^.between the two mtersection points of the normal- 

f etirtr^^""' ^^i""^"^ - determinatTon 6f 

,^film thickness. The contours for these particular 

g^e ambiguity IS seen to be beyond the experimental 
,|rpr, even though the mtersection angles are not large; 
t&?n S,'*'*'™r?;°^ «°"«ideted in the nfx 
Soulin^ "^P^^ Illustrates the usefulness of the 
^£Uque-incidence measurement on a thin fUm. 

^**'"*'0"s fdr n and" 

|Many methods have been developed'.*'." to de- 
«°'^t^'^ts n and ^^from mea^urti 
g^tance^ ancUransmittance T. Because the eqiT 
JiPn3 .7i!(„,&) = R and r(n,A) = T cannot be solved 



analytically even for normal incidence on thin films or 

wTb^rf r ^""^P^^^' S^^Phical mSods 

co^TlH It ^^7^^,."^^*^ decently. Similar methods 
could.be used for this problem, but a large inventory 

wo^H^ ' 7°"/^ }T *° P^^d^'^^d and^usinr?hem 
m^S^H ^""'^ concluded that the bLt 

method is to use a small number of contour graohs^ 

r'^uTv^r"'^ ""'T ^'^^ ?henTo usetl 

comSertime' " ^ 

fn7?Kr''°™^"*^^ inversion technique to be described 
for oblique incidence is similar to ones previously d^ 
If Newton-RapLn Sod" 

of finding the mtersection of contours. A procedure 
that requur^ the explicit calculation of partial 
denvatnres ofR and T with respect to n and AwouM be 
very tedious because of the complexity of theTdaitL 

^hls reqlTm^r^'^"" ""'"^^^^ *^ 

thilL°e^s"i?nf*^hrfif f r ^"'"^ *hat the 
cmcKness d of the film is known by an independent 

SS^t^Slh^^lf ^ ""^^ ^ 4resentThe two 

quantities that have been selected for measurement 
from among the possibilities for R and T a? vS 
angles and polarizations. One starts the procedX^ 
some wave ength with a pair of values ZCST^I 
It^T f^^'^^tf .'>f desired intersection point in 
the plane obtained by inspection of the contours of 
Sec III or by some other means. (In practice one 
TwV^r ""f^ ^ wavelength correspondfng tb a S/X 
£ X^n? """Tr T 3ucceLve steps 

tfon fw ^^V ^'^ preceding sec- 

tion; two intersections may occur in certain regions of 




^'^ff °/ intersection of normal-incidence « and T 
contours with p-polar.zation T contour for 60« incident angle 
d/\ = 0.118. ■ 
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n-k space, the values of uoM should be closer to the 
intersection desired.) One then calculates ^(no,/co) = 
i/'o, 4>(noM) = *o, from the formulas of Sec. II, and 
the partial derivatives at no,fco. Since it is inconvenient 
to calculate these partial derivatives explicitly, they 
will be determined by another method. Assuming 
for the moment that these derivatives are known, one 
next calculates the differences — ^Po, ^ — <^o, where 
^ and $ are the experimental values, and writes 

These two linear equations in the unknowns An,A/c 
can be solved easily. One then writes ni = no + An, 
fci = fco + Afc and inserts these new starting points into 
the procedure. After a small number of iterations, 
A^ and A(t> approach zero and the procedure can be 
terminated with a printout of the corresponding n 
and k values. 

At oblique incidence it is simpler to approximate the 
partial derivatives numerically than to calculate them 
analytically. A least-squares tangent approximation 
can be used, utilizing the matrix method of least- 
squares regression. One starts with the point (no,/co) 
and considers a square of points symmetrically about it 
in the n-fc plane, with corners at (no + Ano, ko + Afco), 
(no + Ano, ko — A/co), (no - Ano, U + A/cq), (no - Anc, 
Ajo — Afco). If Awo and Afco are sufficiently small, say 
of the order 0.01, then the functions ^(n,fc), <^(n,fc) are 
almost planar in this limited region and can be approxi- 
mated by 

^(n,fc) = flo + CLin + ojfc, 



where 



<A(n,fc) ^ bo + bin + h^k. 



where 



ai 



-m^ -{ti- '.-(a 



One next constructs the following matrices: 

1 (no + Ano) (fco + Afco) 



X = 



(no H- Ano) 
(no — Ano) 
(no — Ano) 



(fco — Afco) 
(fco + Afco) 
(fco — Afco) 



^ = 



/ ^(no -h Ano,fco + Afco) \ 
/ ^(rio +• Ano,fco — Afco) I 
I ^(no - Ano,fco + Afco) T 

\ ^(no — . Ano,fco — Afco) / 

' 4>(no + ATiOjfco + Afco) \ 

.(no -h Ano,fco — Afco) \ 

4»ina — Ano,fco 4- Afco) j* 

A (no — Ano,fco — Afco) / 



From least-squares regression theory we have,'*^ de- 
noting the transpose of X by X'', 

a = (X''X)-^X''^, b = (X''X)-^X^4., 



as 




Reading au a,, bi, from the vectors a and b gives the 
partial derivatives. This procedure is repeated at 
each iteration in the Newton-Raphson scheme. Our 
program usually requires three to five iterations for 
convergence of n and fc to one part in 10^ Tested 
against the normal-incidence examples given by Abel^s 
and Th^ye,^^ its results were exactly the same as theirs. 

This procedure may be expanded to include a search 
for the film thickness d if that is not known.^ The 
effect of d on the solution for n and k is illustrated in 
Fig 8. (This figure shows, incidentally, that linear 
interpolation between sets of contours for different 
d/\ is usually satisfactory for a rough estimate of n 
and fc.) The value of d/\ inserted into the two- 
parameter search routine affects the solution for n and 
fc significantly, and if it is too far off, no intersection 
will be found at all in some cases. Measuring d directly 
with accuracy is difficult, especially for thinner films, 
although a standard interferometric measurement gives 
a suitable starting value for the following three- ^ 
parameter search routine. The generalization of the . 
procedure to include a search for d requires measure- . 
ment of three quantities among the various R and T 
combinations, and an estimate of the starting point 
(no,fco,rfo). The generalization", from two to three . 
dimensions is straightforward, and so the details of^;r 
the analysis are omitted. 

Often the measurements are made at many ditferentv 
wavelengths on a single film. The three-parameter fits/ 
at different wavelengths will converge to different values^. 

I- 



n. 
1 




Fig. 8. Loci of the intersections of some particular normal- 
incidence R and T contours for d/\ values from 0.05 to 0.2. 
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of d because of experimental errors in the measure- 
ments, but the physical thickness of the film is actually 
the same at all the wavelengths. Therefore one should 
average the d values obtained at those wavelengths 
where the contour intersections provide accurate solu- 
tions, leaving out those near where tangencies or other 
inaccuracies occur. Then improved values of n and k 
are obtained by carrying out a new two-parameter fit 
using this average value of d together with the best pair 
of R and T values at each wavelength. For example," 
in Ag for light between 0.5 and 6 e V a three-parameter 
fit gave d values which were in the range 310 ± 10 A, 
except near the 4-eV plasma resonance and in the in- 
frared where two of the three contours are nearly 
parallel. The unreliable regions can be recognized 
without reference to the contours by large deviations of 
d from the average or by failure of the three-parameter 
fit to converge at all. The error in the average d value 
can be much smaller than the scatter, perhaps 1 or 2 A. 
This method has the additional advantage of giving 
automatically the (short-range) average thickness over 
the entire area of the film that is exposed to the light 
beam (assuming that the average thickness is constant 
over the illuminated area, which is greater for oblique 
than for normal incidence) . 

V. Comparison of Thin-Film and 
Bulk-Sample Methods 

It is informative to look at some bulk-sample methods 
for determination of n and k in the same terms as our 
discussion of the thin-film possibilities. Such a com- 
parison is necessarily incomplete, since considerations 
of the different instrumental limitations^ are essential 
to an analysis of ultimate accuracy, but it is useful to 
put all these methods into the same perspective. 

Most closely related to the thm-film calculations 
and methods are measurements of R on infinitely thick 
samples at two angles of incidence or two different 
polarizations. 27 Yov this case, T = 0 and i2 = 
|ri2|2. We have already remarked that the 5-polariza- 
tion contours are parallel at all incident angles (except 
near the origin for infinite thickness), so that one can 
rely on a comparison between p polarization at a large 
angle with either s polarization or small incident angle. 
In Fig. 9 are shown the R contours for normal incidence 
and for p polarization at 60° incident angle, for infinite 
thickness. (This is the more useful comparison, be- 
cause the s- and p-polarization contours at 60° become 
parallel near the origin for infinite thickness.) These 
may be compared with the thin-film contours of Fig. 5. 
The bulk contours intersect at small angles for k > 2, 
At 75° incident angle they can be used up to about 
A = 4, except for n < 0.5. Thus we "conclude that a 
measurement of R at two incident angles or polariza- 
tions does not usually encounter the problem of am- 
biguous solutions, but it»Vould yield greater accuracy 
with thinner samples,^ leaving aside problems of sample 
preparation and thickness measurement. These mea- 
surements (especially for k values) are least reliable at 
large k and small n.^^ 



For the Kramers-Kronig analysis, we consider fi^ = 
r exp(t0) in polar form, at normal incidence where fxz = 
{n — 1)/ {n + 1) for both s and p polarization. The 
normal-incidence reflectivity, R = r^, is measured as a 
function of wavelength, and the phase 0 is obtained 
from it by the integral dispersion relation 

^(o,) \nr{oy')d<o' ^ 

. ^ Jo - 
or some variant of it.^s.z^ (The same kind of analysis 
can also be applied to T measurements on thin films.*) 
Then n and k can be expressed explicitly in terms 
of r and 4>: 

n = (1 - r«)/(l + r« - 2r cos^), 

k =^ 2r S\n4>/{1 + — 2r cos0). 

The contours r = const, and ^ = const, are orthogonal 
trajectories of each other, since the function n(fx2) is 
analytic. Thus the problems of tangent contours and 
multiple intersections do not arise when this approach 
is used. In fact the contours are families of circles, 
because the function n{fn) is a linear fractional trans- 
formation. In Fig. 10 the contours for R and 0 are 
shown. The R contours are circles of radius {P^ — 
1)* centered at n = P, A? = 0, where P = (i + R)/ 
0- — R); the 0 contours are circles of radius (Q^ + 
1)* centered at n = 0, A; = Q, where Q = l/tan0. 
From the spacing of the <f> contours, one would conclude 
that the largest uncertainties are introduced in k at 
large k and small n (as in the preceding method). 
Consideration of the accuracy of the integral evalua- 
tion of 0 is beyond the scope of our discussion,, although 
it is important in practice. ^^"^^ 

The Drude method depends on the comparison of 
fia's r, exp(t^s) and fi2^ s exp(i0j,) at oblique inci- 




Fig. 9. Contours of const. R for normal incidence and p polar- 
ization at 60** incidence for infinite thickness. 
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shown in Fig. 12. The angle dependence is not a 
serious problem at 60**. 

VI. Conclusions 

Our general conclusion is that thin-film measure- 
ments at oblique incidence can give accurate values of 
n and k in some regions of the n-k plane where other 
methods are subject to significant uncertainties. Be- 
cause they are especially easy experimentally using a 
standard spectrophotometer in the 1-6-eV region, 
they may well be the method of choice in that wave- 
length range for easily evaporated materials. 



Fig. 10. Contours of const. R and const, phase <^ for normal 
incidence and infinite thickness (Kramers-Kronig method). 



dent angle. One measures ^ and A, where fn^/ fvi = 
(tan^) exp(*A), :or tan = ^j»A* and A = - 
Experimentally ^. and A have usually been measured 
using polarizers and a quarter-wave plate or compensa- 
tor, so that the method has been most used near the 
visible region of the wavelength spectrum. Again n 
and fc can be expressed explicitly in terms of ^, A, and 
^1 (the. angle of incidence). The contours for and 
A are again orthogonal trajectories of each other be- 
cause the relation to n is analytic (although not linear 
fractional). They are shown in Fig. 11 for two dif- 
ferent incident angles. The method must be used with 
large incident angles, since for small angles tan^ 1 
and A TT everywhere in the n-fc plane except near 
n = 1, fc = 0. (In this region smaller incident angles- 
could be useful.) For instrumental reasons a value of 
A near 7r/2 is preferred," subject to some refinements. 
From Fig. 11 we see that satisfying this criterion de- 
mands appreciable changes in incident angle as n and k 
vary with wavelength, but that a relatively large inci- 
dent angle is desirable for most of the n-fc space region 
under consideration. The contour spacing suggests 
that also for this method the errors are greatest at 
large values of fc, but with special care good accuracy 
has been achieved,^^ . 

Finally it should be mentioned that methods m- 
volving p polarization at large incident angle may be 
sensitive to the experimental accuracy of the incident 
angle. These include the bulk reflectivity and Drude 
methods. For example, the dependence on incident 
angle is shown in Fifi 12 for at » thickness, in the 
case n = 1, = 1. If the polarizing angle for this 
n and fc (^52°) is chosen, the problem is minimized, 
but at larger angles it becomes important. For thin 
films, the sensitivity of T with p polarization is also 
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Fig. 11. Contoui-s of const, ratio of amplitudes tan,^ and cons^. ^ 
phase difference Aifor p and s polarization (Dmde method). ^ 
(a) 60° incidence; (b) 75** incidence. - ^ 
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Fig. 12. Dependence on incident angle of R for infinite thickness 
and T for d/\ = 0.1, both with p polarization, in the case n = 

To determine n and k using thin films, the film thick- 
ness must also be determined. It can be obtained 
more accurately from three optical measurements, 
one of which may be p-polarized transmission, than 
by routme mterferometric means. The third mea- 
surement could instead be fl', the normal-incidence 
reflectance from the substrate-metal interface," but 
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the thickness of a thin film even aside from the optical 
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borne metals lie in this region for visible wavelengths. 

Problems of sample preparation have been left out of 
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lems of surface smoothness and oxidation are especially 
important, and can be controlled in evaporated thin 
films. The question of the film structure may be 
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The *V-W' method for absolute reflectance^was first developed by Strong.^ The basic arrangement is shown in 
Figure 1 for (a) sample ' W and (b) reference *.V' measurements. In the sample measurement, two reflections take place on 
the seunple, with an intervening reflection off a mirror. This requires a sample of sufficient size and spatial uniformity. For 
the reference measurement, the sample is removed and the mirror is repositioned by a rotation of 1 80° about the sample front 
surface as shown in Figure 1 (b). Even for opaque samples, careful alignment of the sample and mirror are critical to an 
accurate result. For transmissive samples, the back surface reflections will be displaced according to the angle of incidence 
and the index of the material. The first order contribution from the first back surface reflection is shown by the dashed lines 
and from the second back surface reflection by the dotted line. In the V-W method, the second reflection off the sample will 
thus lead to a larger spread of the reflected light than for absolute methods with a single sample reflection. The use of an 
averaging sphere for detecting the output beam can be used to reduce sensitivity to alignment. 




Figure L Schematic of Strong's V-W method for absolute reflectance measurement The sample measurement geometry is 
shown in (a), in which tvvo reflections take place off the sample and a mirror. For transmissive samples, the first order 
reflected components from the back surfece are indicated as dashed and dotted arrows. For the reference measurement in (b), 
the sample is removed while the mirror flipped and trztnslated. The beam path lengths and angles of incidence ate 
maintained. The absolute reflectance value is obtained by taking the square root of the ratio of sample to reference 
measurements. 

The 'V-N' method is shown in Figure 2. The sample, measured 'N' in (a), does not move for the reference 
measurement. Rather two mirrors are repositioned in the reference 'V* measurement (b). Mirror I is translated without 
rotation and mirror 2 is rotated without translation, so that the beam travels the same path length and is incident on the 
mirrors at the same angles for both measurements. Only one reflection off the sample takes place and the absolute reflectance 
is obtained by a simple ratio of sample to reference measurement result. The additional beams coming from the back surfece 
reflections are displaced from the first surface reflected beam. But the spread is less than in the 'V-W case. As in the case cf 
the * V-W method, accurate aligmnent is critical to high measurement accuracy. Also, an averaging sphere for the detector 
can be used to reduce sensitivity to alignment. ' 
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dashed arrow. For the reference measurement in (b), the sample is removed and the detector is rotated about an axis centered 
on the sample. The absolute reflectance value is the ratio of sample and reference measurements. 

The fourth absolute reflectance method uses an mtegrating sphere as shown in Figure 4 and is detailed in Refea:ence 
1. The sample and reference measurement geometries in (a) and (b), respectively, are shown in top views. The specifics cf 
the locations of the parts of the sphere are important only to the extent that symmetry is preserved between sample^ and 
Inference measurements.' The input and reflected beams, and the sample, reference and entrance port centers, all lie in the 
plane containing the circular cross-section of the sphere shown in the figure. The detector is mounted on the top of the sphere 
and has a field of view (FOV) on the bottom of the sphere. The bafiQe is used to block the sample and reference fi-om the 
FOV, and is located below the plane of the input and reflected beams. The sphere mterior; is coated with a high reflectance 
difluser that acts to uniformly distribute light within the sphere. Besides the reflection off the sample, the only difference 
between the sample and reference measurements is the region of the sphere wall where the beams first hit. These are the 
'specular' and 'reference' regions shown in Figure 4 (a) and (b), respectively, which are symmetrically placed on opposite 
sides of the entrance port. For the infrared sphere discussed in the following section, the ratio of the responsivity of these two 
regions to incident light is very close to 1. 




Figure 4. Schematic of integrating sphere based method for absolute reflectance measurement. Both (a) and (b) are top 
views. The sample measurement is showQ in (a), m which reflection off the sample takes place. The first order back surfece 
reflected component is indicated with a dashed arrow. For the reference measurement in (b), the sphere is rotated twice: first 
about the sample port and second about the sphere center. The baffle, detector and detector FOV are only labeled in (b). The 
" baffle is positioned below the plane of the input and reflected beams. The absolute reflectance value is the ratio of the sample 
and reference measurement. 
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The additional light coming from the back surfece reflections of the sample are disnlacert fr«m fi^f e.,^.» 

umfom coUectmg area of the mtegrating sphere, transmissive samples can be accurately ml^ured In ^hhS ^tw 
attnbutes of the sample such as surface flatness and wedge angle between the front ^dVSs^7^Mc^^J^Tf. ^ 
em>r. m other systems (which are sensitive to aligmnent), have ve,y little effect on r^spSefe ^^S^? ac^S 

A significant difference between the sphere method and the othere described in this section is that th^ h^c 
hemisphencal collection. However, for samples that do not scatter, aU fiibr mXds deliiS Zuld SoS^^^ 
reflectance result. For samples that do scatter, the results will differ depending eS of tht bem"^^^^^^^ ofT 

unplementetions of the other three methods (i.e. solid angles of incidence and measur^e^ More hnnSifdr^^^^ 
fimher m the Section 3. for samples with scatter only the sphere system can provide mISJSgfuVresulte. ^' 

3. QUALITATIVE COMPARISON OF METHODS 

Some of ^e relative advantages and disadvantages of the four methods for absolute reflectance measurement nf tran.™Jcc.„» 
samples becomes apparent in the description of the geometries in Section 2 AroJS>d^^f Jl!^,! ^ 
disadvantages of the four methods is given m Table 1 at le end of this section ^ advantages and 

aHan™.J^° advantages of the integrating sphere method are the insensitivity to alignment and the factois affectine 
ahgmnent. as well as an ability to measure sample scattered light (away from the specular diction) As noted to SeS^ 
factors mcludmg sample thickness, flataess, parallelism, and aligmneni are all important for the non-sXre methods S^d i,' 
result in significant eirors in the detemiination of the sample reflectance. For the sph^ meftod aTof £te fe^S. ™-if ^ 
m small displacements of the light incidem on the sphere wall TTie small iSmente w^niin tn^^^^^ ^'1 
detected signal due to the excellent local spatial unifSmity of gooJ sZ^^e^^^^^^ ^ 

mettoH^*"^^ "P^' K ^ ^T} "'"^^ ^ *° se^itivities tl SiiSent ver oZ a She^ 

method design can be used fen- sorting out sample scattwed light. ^^^rs. nowever, only a spftere 

r^f,^ Became of the low throughput of integrating spheres, especiaUy in the infrared, the measurement of verv low 
mftSf P '^1* T***'* ^^''"^ ^ » signal-to-noise ^tio, as compared To Ae ^"n^ ^d gon'omSe^ 

T^^JZ^:^^'^'"^^?' viability of the integratmg sphere method is dependent on Te use T FT 
t^S^ £ ^'^ sigaal-to-noise advantages of FT instruments, the noise component of the mncertalL woSd li 

of the ahiit°?f " ^^-"^ ^^^^^'^ '^P'^ scattering light to some degree in some spectral region Verification 

uL^ to n.^ h^ff 't^ important as accurately measming tiie scatter tiiat may exist.» The integ^K sphSr^ S 
used to perform both fimctions. This can be done for bofli reflectance and tiansmittence For the meSmfnt of cL^r 
reflectance, a com|,ensating wedge can be used so tiiat tiie specularly reflected beam to Fig^ 4 (a) eXSTphl^^^^ 
?o^»f^'* °f scatter to transmittance. tiie spheri can be fotLd abo^th^siSpr^ £ 

nomalttansmittance position so that the specularly transmitted beam exits &e sphere through °" 

Because of the pair of reflections off the sample to the 'V-W method, high reflectance samnle<i can in nrm^i.,i» 
measured more accurately than with other methods. When the square root ohhe rerulfj^erfte c^i^nonH^/^^^^ 



reflectance obtained from handbook values of the index of refraction, 
associated with the measurement made at 4 cm" resolution. 



The local structure in both spectra is the result of noise 
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fliigure 5. Absolute reflectance and transmittance j of silicbn window sample, as measin-ed.by the infrared integrating sphere 
[system, and compared with calculated values using handbook n and k valu^ 

The absorptance determined from the reflectance and transmittance experiniental data is plotted in Figxu-e 6. For the 
fispectral range of measurement, lio absorptance. should be seen. Hence the deviation .froni 0 absorptance Js .§L gauge..ofjg?e 
ineasurement error of the reflectance^and transmittal level of ±0.002 in Figure 6 (ignoring the noise) is^, a reasonable 

estimate of the uncertainly of both the reflectance and transmittance measurements. Based only on these results this 
conclusion might be unwarranted, sincj? conceiyaibly, errors m the two measurements, could.^be of opposite sign and tend to 
cancel. However, numerous IR materials have been examiiied using the sphere system with reflectance values ranguig up to 
0.998 and down to 0.05 and corresponding transmittance values down to 0,002 and up^ to .0.95. The results of all of these 
measurements have been consistent with a maximum cuniulative error in the range of ±0.003. 

^. Thus the transmittance measurement result can be used to indirectly evaluate the accuracy of the reflectance result for 
tr^missive samples. In principle designs incorporating transmittance capability are feasible for all four methods. Although 
not commonly found on 'V-W accessories, transmittance nieasurement in the same setup is posisible. The sample 
tfansniittance measurement can be made by translating the sample from its position in Figure 1 (a) so that the sample's edge 
is lined up vertically with the center of the muror and the mkror rnoved to its position m Figures! (b). Then the input beam 
can pass through the sample and the mirror reflected beam in Figure 1 (b) can clear the sample's edge. 

\<^' The 'V-N' geometry shown in Figure 2 does not accommodate a modification for transmittance ine^urement 
V However a hybrid design which incorporates features of the *V-W' method can. The hybrid design sample measurement 
consists of that shown m Figure 2 (a), except with the sample and mirror 1 interchanged. The reference measurement is 
performed with mirror 2 rotated as in Figure 2 (b), but with mirror 1 relocated and the sample renioved as in the ' V-W case. 
" In. the hybrid reference measurement geometry, the sample transmittance can be measured by leaving the sample in position. 

vf For the gonion^ter method, the transmittance measurement geometry is achieved by translating the sample back 
^ into place in the referenc? geometry shown in Figure 3 (b). For the sphere method, dib transmittance measurement geometry 
; is'achieved by a pau* of rotations placing tiie sample in the beam path of Figure 4 (b). This results in beam inciderice on the 

same specular region as in the reflectance case. 
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4. DIRECT COMPARISON OF GONIOMETER AND SPHERE METHODS 

Two methods for absolute reflectance described in this paper have been implemented on a FT spectrophotometer for near and 
mid IR operation (1 jxm to 19 }im).' These are the integrating sphere and goniometer methods described in Sections 2 and 
3. They are shown in photograph form in Figures 7 and 8. 

The integrating sphere system shown in Figure 7 has been designed and constructed for absolute reflectance 
measurement of both diffuse and specular samples. The sphere is 15 cm in diameter with an entrance port of 33 mm diameter 
and sample and reference ports of 22 mm diameter. Baffles inside the sphere shade the sample and referraioe ports from the 
detector port on the top, and the detector FOV on the bottom, of the sphere. The inside wall of the sphere is coated with a 
nearly Lambertian diffuser that has a directional hemispherical (diffiise) reflectance of 0.9 for the infrared spectral range (plasma 
sprayed metal over-coated with gold). The IR beam from the FT Spectrophotometer is focused onto the sample at an 
mcident angle of 8*" and an f/5 cone. Spot sizes on the sample can be adjusted as required. 




Figure 7. Rear view of infrared integrating sphere system for absolute reflectance and transmittance. The sample and refaence 
"ports are framed above and below by adjustable holders designed for use with transmissive samples. The dewar of the MCT 
'detector can be seen at the top of the sphere. A pair of inotorized rotation stages is located underneath the sphere for orienting 
it to the various positions for transmittance, reflectance, reference and other measurements. 

An Hg:Cd:Te detector with nonimaging concentrator optics is mounted on the top of the sphere shown in 
Figure 7.*^ The entrance port for the input beam in reflectance mode is located on the back of the sphere in this photograph. 
The sample and reference ports can be seen in the foreground, with adjacent adjustable holders. These are designed to avoid 
^ obstruction to the sample transmitted light in the reflectance mode and to the input beam in the transmittance mode. Below 
' the sphere are two rotation stages. The bottom stage has its axis of rotation centered on the input beam focus. The top stage 
' holds the mtegratmg sphere centered on its rotation axis and is mounted on the bottom stage off center. The top stage is 
' used to select which sphere port is placed at the beam focus, and the bottom stage is used to vary the incidence angle of the 
•* beam on the selected port (from 0° to 360°). Hence reflected light can be measured in the orientation shown m Figure 7, 
' while transmitted light canf be measured with a 154^ rotation of the bottom stage. Details of the sphere design and 
' ' perfomiance beyond those presented here can be found in References 1, 7 and 12. 
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With the addition of a purge enclosure for both systems, made prior to the sapphire comparison measurements, these 
^ uncertainty levels apply to the entire spectral region. 
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Figure 9. Direct comparison of integrating sphere and goniometer absolute reflectance measurements of a protected Ag mirror 
in the infrared spectral region. Structures from 2.6 |xm to 3.0 ^m, at 4,2 urn and from 5.6 iLim to 7.0 |im are due to 
atmospheric absorption effects not related to the Ag mirror. 

A sapphire infared window was also characterized using both the integrating sphere and goniometer systems. Both 
the reflectance and transmittance were characterized on both systems. The results are shown in Figure 10. A purge enclosure 
was in operation for these measurements. The sajiphire window has an anti-reflection coating optimized for 4 \im 
wavelength. Due to the coating, the transmittance and reflectance in Figure 10 vary considerably from uncoated sapphire 
values. The agreement between tiie sphere and goniometer measurement results appears good. A closer look at the results is 
. shown in Figure 11. * 

For each system the absorptance is calculated from the transmittance and 'reflectance data shown in Figure 10. The 
absorptance values nearly span the complete range from 0 to 1. To compare the two method's results, the difference between 
the two absorptance spectra is plotted m Figure 11. Two primary features stand out. The feature at 8 |JLm is not understood 
at present, but is under study. The sinusoidal variation for shorter wavelengths is attributed in part to the angular 
dependence of the anti-reflection coating. Samples with a strong sensitivity to angle of iricidence will naturally have larger 
uncertainties associated with their measurement Both the difference in the beam geometries of the two systems, and the 
slight variation in angle of incidence for the transmittance and reflectance sample measurements of the two systems can result 
in larger 'errors' in the calculation of absorptance. More accurate angle of incidence aligimient in both vertical and horizontal 
planes, combined with corrections for the difference in input beam geometries is expected to reduce the size of the variations 
seen in Figure 11. 
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5. CONCLUSIONS 



Four methods for the absolute measurement of reflectance have been described and evaluated in terms of the characterization of 
transmissive materials. Certainly these materials present more difficulty for accurate measurement than opaque materials do. 
^ A qualitative coniparison of the four methods has been presented along with a quantitative comparison of infiared 
implementations of two of the methods: the integrating sphere and goniometer methods. In most cases alignment errors 
dominate the measurement uncertainty. Where feasible, integrating spheres can be used to significantly reduce these errors, as 
■ in the case of the integrating sphere method. The incorporation of transmittance measurement under the same beam geometry 
conditions is an important tool in the evaluation of the reflectance results and measurement uncertainties. 

Two methods for absolute reflectance described in this paper have been implemented on a FT spectrophotometer fir 
the infrared spectral range. In the direct comparison of these two methods using two types of samples, the agreement between 
the infrared integrating sphere and goniometer systems is reasonably good and within the estimated levels of uncertainty. 
^ The planned addition of larger area, more spatially unifonn detectors or averagmg spheres should result in improved 
performance of the goniometer system. 
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